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By Murray Tobak 
Ames Research Center 


SUMMARY 


A method is presented for deducing certain nonlinear differential 
equations from their solutions. The method,, relevant to equations of the type 
x + €[xg(x 2 ) + xk(x 2 )] + v 2 x = 0, enables the deduction of explicit forms for 
the nonlinear elements g(x 2 ) and k(x 2 ) . It should apply whenever the solu- 
tion can be assumed to be describable by the first approximation of the 
Krylof f-Bogoliuboff method . 


INTRODUCTION 


Many fields of mechanics require the study of oscillatory systems 
governed by nonlinear second-order differential equations. A wide class of 
such equations yields to treatment by the method of Kryloff and Bogoliuboff 
(ref. l) . In particular , the first approximation affords a ready means of 
finding a solution capable of revealing the main features of the oscillatory 
system . 

It does not seem to have been appreciated that the first approximation 
also lends itself readily to a problem which is the inverse of the usual one; 
namely, given the solution, find the governing differential equation. This 
problem arises continually, for example, in experimental studies of aircraft 
dynamics. Here, the result of an experiment is often a time history of an 
oscillation. The analyst must extract from this history the form and magni- 
tude of the aerodynamic quantities which are judged to influence the motion. 

In effect, therefore, the analyst must reconstruct the differential equation 
which has yielded the observed record. The purpose of this note is to show 
how the first approximation of the Krylof f-Bogoliuboff method may be adapted 
to this end for a certain class of equations . For equations of the type 
X + e[xg(x 2 ) + xk(x 2 )] + V 2 x = 0, it is shown that, given the solution, 
explicit forms for the nonlinear elements g(x 2 ) and k(x 2 ) are deducible from 
inversions of Abel integral equations. 

A note on precedents: The method proposed here is of such simplicity 

that it raises the question whether it has not been proposed before . A search 
for precedents has disclosed one, applicable, however, to another class of 
equations. In reference 2, G. Plato has shown by a heuristic argument (con- 
firmed by the method proposed here) that the form of the nonlinear element 
f(x) of the equation x + ef (x) + v 2 x = 0 also can be deduced from an Abel 
integral equation. It is easy to envisage other classes of equations which 
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will yield to the same treatment. It is probable that similar methods, appli- 
cable to still other classes of equations, have been or remain to be 
discovered. 


ANALYSIS 


The First Approximation 

The analysis of reference 1 is concerned with differential equations of 
the form 


x + v 2 x + ef(x,x) = 0 


(i) 


where € must be a small quantity. Since, with e sufficiently small, x 
will not deviate greatly from the harmonic solution x = a sin(vt + 9) , a 
solution is sought in the form 

x = a(t)sin ^(t) (2) 

where 

t(t) = vt + <p(t) ( 3 ) 


It is found that, to a first order in e, a and ijr must satisfy the relations 


e r 23t 

a = / f(a sin ty, av cos cp)cos cp < 3.9 

2 itv J 


(h) 


\jr = v + 


2jtav 


L 


2Jt 


f(a sin cp, av cos cp) sin cp dcp 


(5) 


When equation (l) is given explicitly and the solution for x is desired, 
equation (2) provides a solution, a(t) and \|r(t) being determinable directly 
from equations (4) and ( 5 ) • 
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The Inverse Problem 

Now consider the inverse problem. Let it be assumed that the solution is 
known in the form of equation ( 2 ); that is, a(t) and ^(t) are known functions. 
It is desired to find the form of equation (l) under the assumption that the 
solution is adequately described by the first approximation. This assumption 
enables one to seek a solution for f(x,x) through equations (4) and ( 5 ). 

Consider a class of problems for which it is possible to specify the form 
of f(x,x) to the extent shown in equation ( 6 ) 


f(x,x) = xg(x 2 ) + xk(x 2 ) 


(6) 


This form is chosen principally on physical grounds, anticipating in particu- 
lar its pertinence to aerodynamic applications. In reference 3 it is shown 
that, at least for low reduced frequencies, equation ( 6 ) is a general form for 
the aerodynamic pitching moment when x is identified with angle of attack. 
•The specification that g and k be even functions of x again envisages 
aerodynamic applications. 

Inserting equation (6) in equations (4) and ( 5 ) gives 


e r 27t 

a = - / cos <pt(av cos $)g(a 2 sin 2 cp) + (a sin cp)k(a 2 sin 2 Cp)]dcp 

o (7) 


\|/ s= v + 




2rtav 


sin cp[(av cos cp)g(a 2 sin 2 cp) + (a sin cp)k(a 2 sin 2 q>) ]d<p 


( 8 ) 


In equation (7)> the term involving k, being odd about cp = jt, contributes 
nothing to the integral and may be dropped. In equation ( 8 ), on the other 
hand, the term involving g may be dropped for the same reason. After the 
change in variable £ = a sin cp, equations ( 7 ) and ( 8 ) take the form 


dt 


(a 2 ) = - 



i 2 g(i 2 )di 


( 9 ) 
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( 10 ) 


ijr = v + 


2s r a 

7a 2 - I 2 


Ttva" 




With a and \|r being known functions, equations ( 9 ) and (lO) are a pair of 
integral equations, the first involving the unknown function g(£ 2 ) alone, 
the second involving the unknown function k( | 2 ) . 


Inversions 


Consider first equation ( 9 ) • Put a 2 = 0, £ 2 = A, and let 


g(A) 

a/A 


= 9(A) 


(11) 


so that 


e = - 



o 


( 12 ) 


Since 0 is known, it may be assumed that 0 is also known, and moreover, 
that 0 can be represented, albeit often only numerically, as a function . 
of 0. 1 Writing 


0 = P(0) 


(13) 


and differentiating with respect to 0 in equation (l2) yields 



0 q(A)dA 
\l 0 ~ A 


(i*0 


Equation (l4) is recognized as Abel's integral equation, which has a particu- 
larly simple inversion (see, e.g., ref. 4). It is 


1 In practice, it may be necessary to have a number of records, covering a 
range of initial conditions, to ensure an adequate representation of this 

dependency. C han ges in initial conditions should only extend the range of the 
dependency without destroying its uniqueness; the function 0 = F(0), in other 
words, should be independent of the initial conditions. 
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(15) 


cl(0) = - 


A f 

&9J 


F* (A)dA 
n fe - A 


This determines q, and hence g, through equation (ll) . 

Equation (lO) is treated similarly. Again putting a 2 =9, £ 2 = A, and 
letting 

JK k(A) = p(A) ( 16 ) 


yields 


■4 


(17) 


With 


ve(^ - v) = g(q) 


( 18 ) 


the inversion of equation ( 17 ) is 


= l A. I " 9 Q ( A ) dA 
p( } ~ € MJ 


(19) 


This determines p, and hence k, through equation (l 6 ) 


Example 

To illustrate the method, we adopt a problem considered in reference 3- 
Let us assume that from one or several records of the oscillatory system, it 
has been found possible to represent the solution x(t) as 


x = 


■fl - e- £t (l - VO 


sin 


jt - 1 


|i In 


1 + 


x r 


(e tt - 1) 


+ ♦, 


>} 


( 20 ) 


Note that x approaches a limit motion as t -*■ oo. From equation ( 2 ) we 
identify a, and v as 


■i. 
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e 


4 


a 


2 


1 - e" et (l - 4/x 0 2 ) 


'S 


t = t 


3 

| H irx 





( 21 ) 


v = 1 

J 

In the determination of g ; it is required first that 0 "be written as a 
function of 0. A straightforward calculation yields 


0 = F(0) = - | 0(0 - ^) 


( 22 ) 


One notes that, as required, F(0) is independent of the initial conditions. 2 
Inserting F’(0) in equation (15) and performing the indicated operations 
gives 


q.(e) = j=i° - 1) 


(23) 


whence, from equation (ll) 


g(0) =0-1 


(24) 


The determination of k proceeds similarly. Writing \|r - v as a function 
of 0 yields 

G(0) = - | P£0 2 (25) 

Inserting G in equation (19) leads to 

vie) = -p 0 3/2 (26) 


^The advisability of having a number of records, covering a range of 
initial conditions, is emphasized in the case of limit motions. Here, for 
certain initial conditions, 0 may be constant over all or a major part of the 
record. In the present case, for example, when x Q = 2, 0=4. If this were 
the only record available, one might erroneously conclude that 0, and hence 
g, were identically zero for all 0. 
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so that, from equation (l6) 


k(0) = -ia© 


Hence, in equation (6) 


f(x,x) = x(x 2 - l) - ^x 3 

so that the equation of motion is, correctly (e-g., ref. 3) 


x - ex(l - x 2 ) + x(l - epx 2 ) = 0 


( 27 ) 

(28) 


( 29 ) 


DISCUSSION 


The validity of the proposed method for determining the form of the 
governing differential equation from its solution rests entirely on the 
assumption that the solution be describable by the first approximation of the 
Kryloff-Bogoliuboff method. Unfortunately, it is not invariably the case that 
the appearance of the solution itself is sufficient evidence to justify the 
assumption a priori; an a posteriori justification is necessary. 

Having obtained a form for the differential equation on the basis of the 
assumption, one should determine that over the range of amplitudes and fre- 
quencies of the solution from which the equation has been derived, the 
inequalities 


max | eg(x 2 ) | « v 


max | €k(x 2 ) j « v 2 


( 30 ) 


have been satisfied. Failure to satisfy the inequalities must be taken as an 
indication that the derived equation may be incorrect. It is advised that in 
this case, a new evaluation of the nonlinear elements be made on the basis of 
a smaller range of amplitudes and frequencies of the solution, such that the 
inequalities are satisfied. If analytic expressions are matched to the non- 
linear elements so determined, the resulting differential equation then may 
be integrated numerically over an extended range of amplitudes and frequencies 
to check whether the solution it yields matches that part of the solution 
which had to be excluded in determining the equation. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 9 , 1965 
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